Frequency- and transverse wave-vector-dependent spin Hall conductivity in 
two-dimensional electron gas with disorder 
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We determine wave number q and frequency u) dependent spin Hall conductivity o~y x (q, ui) for a 
disordered two dimensional electron system with Rashba spin orbit interaction when q is transverse 
to the electric field. Both the conventional definition of spin current and its new definition which 
takes care of the conservation of spins, have been considered. The spin Hall conductivitivities for 
both of these definitions are qualitatively similar. o-y X (q,ui) is zero at q — 0, uj — and is maximum 
at q = and at small but finite uj whose value depends on different parameters of the system. 
Interestingly for uj — > 0, o~y X (q) resonates when A ~ L ao which are the wavelength (A = ^ir/q) of 
the electric field's spatial variation and the length for one cycle of spin precession respectively. The 
sign of the out-of-plane component of the electrons' spin flips when the sign of electric field changes 
due to its spatial variation along transverse direction. It changes the mode of spin precession from 
clockwise to anti-clockwise or vice versa and consequently a finite spin Hall current flows in the bulk 
of the system. 

PACS numbers: 72.25.-b, 71.70.Ej, 72.15.Lh 
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I. INTRODUCTION 

One of the primary goals today in spin based 
electronics 1 is the generation of spin current. Recent 
realization of spin-Hall effect (SHE) 2-4 in semiconduc- 
tor systems is certainly a very significant achievement in 
this direction. This is a phenomenon for electrical gen- 
eration of spin: a charge current along its transverse di- 
rection induces a spin current whose polarization is per- 
pendicular to the plane formed by these two currents. 
This phenomenon was predicted 5-7 long back due to the 
spin-asymmetric 'skew-scattering' mechanism and spin 
dependent 'side-jump' mechanism which are collectively 
called extrinsic mechanism because the spin-orbit inter- 
action (SOI) is disordered in this case. This mechanism 
is also responsible for anomalous Hall effect (AHE) 8 ' 9 
in ferromagnets. However uniform (pure) SOI which 
is intrinsic 10 , also causes AHE. Similar intrinsic mech- 
anism due to the SOI in hole-doped semiconductors 11 
and two dimensional electron gas 12 in semiconductor het- 
erostructures have been predicted to give rise to nondis- 
sipative spin Hall conductivity (SHC). 

The spin accumulation observed in n-doped GaAs 2 is 
believed to be extrinsic in origin because of small spin ac- 
cumulation and its directional independence on the elec- 
tric field. On the other hand spin accumulation in two 
dimensional hole gas (2DHG) is large 3 and hence is sug- 
gested to be intrinsic in origin. These experiments do not 
measure the spin voltage or the spin current, however the 
technique developed in observing charge accumulation 13 
at the transverse edge due to spin current which is called 
inverse SHE, could be useful to measure spin Hall cur- 
rent. Nevertheless estimated SHC from the observation 
of spin accumulation 2 in n-doped GaAs is in good agree- 
ment with the calculated 14,15 SHC in extrinsic mech- 
anisms. An effective two band cubic (in momentum) 
Rashba model describes the electronic states in 2DHG 



well. Even in presence of disorder, the SHC in 2DHG 
has nonzero 16 intrinsic contribution. 

The responsible mechanism for SHE in two dimen- 
sional electron gas (2DEG) which we study here is par- 
ticularly not clear yet. The electronic states in the 2DEG 
formed in semiconductor heterostructures can be well de- 
scribed by the Hamiltonian 



H = — <t + A(fc y (7i - k x a 2 ) 
2m 



(1) 



where k and m are the momentum and mass of the 
electrons respectively, A is the Rashba spin orbit cou- 
pling strength 17 , cto is the unit matrix, and <7j are the 
Pauli matrices. (We have set the unit % = c = 1.) 
Sinova et al 12 predicted universal SHC o s yx — e/Sir for 
such systems using conventional definition of spin cur- 
rent j Q = i {s a , i>k} where group velocity = Vki?o 
and s a — <7 a /2 is the ev-th (a = 1 — 3) component of 
spin. However after a prolonged debate 18-28 , the con- 
sensus arising from various methods of calculations is 
that <jy X = in presence of disorder, no matter what 
its strength is. These studies include calculation of ver- 
tex correction 19 in Born approximation, using Keldysh 
formalism 20 ' 21 ' 27 for any value of lifetime r, using Kubo 
formula analytically 22 ' 23 ' 28 and numerically 25 ' 26 , and 
Boltzmann transport equation approach 24 . The equation 
of motion for spin projected on the plane is dt{<Ji,<T2) = 
— 4mA(j 3 , j 3 ). A very unique feature 18 ' 22 of the linear 
Rashba model is that d t <J2 is proportional to j 3 . It sug- 
gests zero spin Hall current in the steady state 22 . This 
simple argument describes the vanishing a s yx for such sys- 
tems. A similar argument also describes spin-spin-Hall 
current 28 since j 1 = — j 2 for steady state derivable from 

the equation dt<Jj, = 4mA(j' 1 + j 2 ). All these results 
are obtained from the above conventional definition of 
j a which is not conserved. The new definition of con- 
served spin current proposed by Shi et al. 31 gives rise to 
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vanishing total a yx for short-ranged 5(r) impurity poten- 
tial and for long-ranged potential upto first order Born 
approximation 27 . 

Is then the 'intrinsic' mechanism really absent for spin 
Hall effect in 2DEG? In a disordered 2DEG, an in-plane 
applied magnetic field may lead to the nonvanishing in- 
trinsic SHC 32 due to Zeeman coupling. Further the inter- 
play of Zeeman coupling with different spin-orbit inter- 
actions may also lead to finite SHC 33 in a pure system. 
In this paper, we calculate SHC using Kubo formula at 
finite frequency to and momentum q transverse to the 
applied electric field within the intrinsic mechanism in 
a disordered 2DEG with no applied magnetic field. We 
find that even in the static limit, SHC is nonvanishing 
and hence the presence of 'intrinsic' mechanism for spin 
Hall effect in 2DEG is demonstrated. 

The paper is organized as follows. In the next section, 
we calculate frequency and transverse momentum depen- 
dent SHC in a disordered 2DEG with Rashba SOI using 
Kubo formula with the conventional definition of spin- 
current. The contribution of spin torque to the SHC 
is also calculated and this contribution, shown in Sec- 
tion III, is qualitatively similar. We find that SHC is 
reasonably high at some range of frequencies and mo- 
menta. Particulary interesting case is for static but spa- 
tially varying electric field: The SHC resonates when the 
wavelength of the spatial variation of the electric field 
matches with the spin precession length. A simple mech- 
anism for this "anomalous" spin Hall current in 2DEG 
is described in Section IV. Section V is devoted for an 
experimental proposal to test this novel mechanism, dis- 
cussion and summary. 



II. SPIN HALL CONDUCTIVITY 

The spin Hall current for an electric field E(q,w) at 
the wave vector q transverse to the direction of E and at 
the frequency u, j y (q,oj) = <r° x (q,uj)E x (q,uj). The spin 
Hall conductivity is nonlocal; spin current at a position 
r depends on the electric field surrounding it: j^(r, u>) — 
J dr'(jy X (\r — r'|, uj)E x (r' , u>). Using Kubo formula 29 , we 
find the transverse spin Hall conductivity 



{i°(k+|) + jO(q,c)}^ +q H] (2) 



with 



Here retarded (advanced) Greens function for an energy 
e can be written as 
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FIG. 1: (Color online) Transverse spin Hall conductivity ay X 
in the unit of 8 2 (e/2iv) as a function of q/(2m\) (horizontal 
axis) and ui/eF (vertical axis). The parameters A = 10, 8 — 
0.4 (upper- left panel), A = 20, 8 = 0.4 (upper-right panel), 
A = 10, 8 = 0.8 (lower-left panel), and A = 20, 8 = 0.8 
(lower-right panel) are considered. 



Equations (2) and (3) together describe sum over infinite 
series of ladder diagrams. We solve the matrix equa- 
tion (3) numerically and then using Eq. (2) we calculate 
&y X (q, ui) when E \\ e x and q || e y , i.e., transverse <Jy X (q)- 
For a system with Fermi energy cf and spin-splitting 
energy 2Akp with kp being the Fermi momentum, we 
choose two parameters A = tpT and S = 2\kpT compar- 
ing with disorder broadening 1/t. 

We show (Jy X {q, uj) for A = 10 and 20, and 5 = 0.4 and 
0.8 in figure 1. The standard resonances occur at finite u> 
and at zero or very low value of q. The maximium value 
of Uy X (q, ui) is almost proportional to S 2 and it decreases 
with the increase of A. One common interesting feature 



for different combination of the parameters to notice is 
that the value of <jy X (q,u>) is not small for u> — > and 

A = 10 



Figure 2 shows a yx 
and S — 0.1, 0.4, 0.8 and 1.2. These choices of S corre 
where L 
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FIG. 2: (Color online) Spin Hall conductivity crjL (dot-dashed 
and dashed lines) for the conventional definition of the spin 
current vs. q/(2m\) for 8 = 0.1,0.4,0.8 and 1.2 from top to 
bottom and for a fixed value of A = 10. af, T is in the unit of 



(Ty X is indeed independent of A as we have checked 
The maximum value of a yx occurs 



ed 2 /2n. 

for an wide range of A > 1 

at q/(2mX) ~ 1, The width of the peak in a yx is larger for 
smaller values of 5 and the value of the peak is larger for larger 
values of S. Solid line represents the analytical expression (9) 
of a yx for small S. 



traversed by an electron while its spin precesses by one 
cycle. We have checked that a yx (q, 0) is independent of 
A for an wide range of A > 1 while S is fixed and is al- 
most proportional to S 2 . a yx is zero at q = as we know 
from various calculations 19-28 , and then it gradually in- 
creases with q and form a peak around q ~ 2mA before 
it vanishes asymptotically. The position of the peak is 
almost independent of 6 but does depend on A. Since 
q/(2mX) = L so /A, the spin Hall current is maximum 
when A ~ L so . In the limit of small disorder broadening, 
i.e., for large 5, <J yx {q) resonates exactly at q = 2mX as 
in the case of 5 = 1.2. Choosing different values of <5 for 
a fixed value of A implies different values of L so . Larger 
the value of S means smaller L so . As the value of L so 
becomes smaller, the decrease in SHC will be faster from 
its peak value for both increase and decrease of A. This is 
the reason for narrower width of the SHC peak for larger 
values of S as we see in Fig. 2. The peak value of <J yx (q) 
is larger for larger 5, i.e. for larger A as well as r. 

To demonstrate the resonance in cr yx (0, q) analytically, 
we calculate J x in Eq. (3) for q x = and q y = q. We 
sum over infinite series of ladder diagrams starting with 
the contribution from ladder with just one bar, 



only Ji and J 3 survive at q x = and they are 
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Since the resonance occurs in cr yx (0, q) at q ~ 2mA <C hp, 
we may wish to evaluate lap up to quadratic in q and 
hence the relevant components are J22 = 1 — (A/m)rg 2 — 
S 2 /2, I 33 = 1 - (A/m)rg 2 - S 2 , I 32 = -I 23 = 2i5Aq/k F 
for 5 1. In this approximation, J 9 = eX8 2 /2 and 
J 3 = —iqeS 2 /4m. Using Eq. (2), we thus find 



<x(0,?) 



(e/27r)<5 2 g 2 
(g 2 + l)(<f + 2) 



1 



(q 2 + 1) 



(9) 



where q = q/(2mX). As we have seen in our numerical 
evaluation, a yx is independent of A and proportional to 
S 2 for small S. The expression of a yx (9) is graphically 
shown in Fig. 2. It agrees well with the numerical evalu- 
ation at low q. The discrepancy at higher q is expected 
as we have evaluated I a p, J® and J3 analytically up to 
quadratic in q only. Nevertheless the analytical expres- 
sion (9) explicitly shows the resonance in a yx (0, q). 



III. SPIN HALL CURRENT FOR SPIN TORQUE 

The above calculation of a yx has been performed using 
the conventional definition of the spin current. We now 
consider the new definition of the spin current which 
is defined to satisfy continuity equation d t S 3 + 'V ■ J 3 = 
and can be expressed as the sum of J 3 and spin torque 
dipole density P T , i.e., J z = J 3 + P T as proposed by 
Shi et aP 1 . Here S 3 and J 3 are the spin density and 
conventional spin current density operators respectively. 
Further the spin torque density operator is expressed as 
r 3 (r) = — V • P T (r) since the average torque density van- 
ishes in the bulk of the system. The second quantized 
form of the spin torque is 



r 3 (q,<) = £ Cl a (t)T af3 (k+ 3)C k+qi/3 (i) (10) 



k,a,/3 



where (t) is the electronic creation operator of mo- 
mentum k and spin a (up or down) at time t, and the 
spin torque operator f (k + 3) = A(k + 3) • <r. We define 
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Expressing J x = Y^, a =o ^a a a and summing over geomet- 
rical series obtained from ladder diagrams we find that 
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FIG. 3: (Color online) Total transverse spin Hall conduc- 
tivity Gy% in the unit of 5 2 (e/2iv) as a function of q/(2m\) 
(horizontal axis) and cu/eF (vertical axis). The parame- 
ters A = 10, 5 = 0.4 (upper- left panel), A = 20, 6 = 0.4 
(upper-right panel), A = 10, 8 = 0.8 (lower-left panel), and 
A = 20, S = 0.8 (lower-right panel) are considered. 



such that T3(0, q, ui) = Q^.°(0,q,u>)E x . Therefore the ex- 
tra part of transverse spin Hall conductivity for an ap- 
plication of electric field along ^-direction is 



9fQJ°(g,w) 



(12) 



The total (conserved) transverse SHC is then <jy S x {q,u)) — 

a yx (q,u))+ OyJ (q,u)). The values of o c y s x (q, uS) for A = 10 
and 20, and 5 = 0.4 and 0.8 are shown in Fig. 3. The 
total transverse SHC is qualitatively similar to the same 
calculated for conventional spin current shown in Fig. 1. 
Figure 4 shows &y X (q) at zero frequency for A = 10 and 
6 = 0.1, 0.4, 0.8, and 1.2, although the value of A is 
immaterial. The only relevant parameter is 5. The peak 
in SHC occurs at q/(2mX) ~ 1 and the width of the 
peak is larger for smaller value of S as in the previous 
case shown in Fig. 2. Also the maximum value of &y S x (q) 
is larger for larger 5 and it is almost two times that of 
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FIG. 4: (Color online) Total spin Hall conductivity er^(g) as 
a function of q/(2m\). Other quantities and descriptions are 
same as in Fig. 2. 
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FIG. 5: (Color online) Electric held along a>axis and its varia- 
tion along y-axis with wavelength A. Electrons moving along 
y-axis will have spin precession in the y — z plane. L so is 
the spin precession length which is the length traversed by an 
electron while its spin precesses by an angle 2n. Lines with 
arrow indicate the direction of spin while it precesses. The 
mode of rotation of spin changes when the sign of the electric 
field changes. Therefore the spatial sign change of the electric 
field induces a net out-of-plane spin current in the transverse 
direction. The spin current will be maximum when A ~ L so . 



IV. MECHANISM FOR ANOMALOUS SPIN 
HALL CURRENT 

The time derivative of the charge current operator is 
jiven by 



ftj£= -5 t k-4emA 2 (e z x j£) 
m 



(13) 
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Therefore in the presence of external electric field E, the 
steady state equation becomes 

e l E = M + Aem\ 2 {e z x <j 3 }) (14) 
m t 

where e z is the unit vector along z-direction and angu- 
lar brackets represent average value of the quantity in- 
side angular brackets. Therefore the electric field creates 
not only the charge current but also the spin current for 
an electron with momentum k. Although some electrons 
move perpendicular to the electric field, their spin precess 
about Ilk = (k y e x — k x e y ) / \ while moving and hence no 
net spin current flows in the bulk of the system. This is 
an alternative description for vanishing a yx in the system 
for static and uniform electric field. However the situa- 
tion alters when the electric field is static but nonuniform 
along its transverse direction. The basic physics behind 
the anomalous behavior for spatial variation of the elec- 
tric field is described in Fig. 5. If the electric field is 
applied along x-direction, due to the wave propagation 
along y-direction, sign of the electric field changes alter- 
nately along y-direction with the wavelength A = 27r/q. 
The spin of an electron moving along y-direction will pre- 
cess in the y—z plane with the precession length L so . 
However at the position where the sign of the electric 
field changes, the clockwise (anticlockwise) precession 
will change into anticlockwise (clockwise) precession as 
in Fig. 5. This is because the sign of the z-component of 
the spin changes as described by Eq. (14). The mode of 
spin precession for the electrons moving along negative 
y-direction will be exactly opposite and therefore there 
will be a net z-polarized spin-current along y-direction. 
The spin current will be maximum when L so ~ A and 
it will sharply fall for the change in A either way. This 
argument for anomalous spin hall current is somewhat 
similar to Pippard's 30 ineffectiveness concept for anoma- 
lous skin effect in metals but the fundamental difference 
is that the former occurs for A ~ L so while the latter is 
due to A <C I (£ is the mean free path of an electron). 
The present picture is ofcourse valid for I < L so . 

We have determined the transverse wave number q de- 
pendent spin Hall conductivity by using Kubo formula in 
the previous two sections and show here how the basic 
physics presented in Fig. 5 describes the anomalous spin 
Hall current. On the other hand, there will be no longitu- 
dinal spin Hall current (when q || E) for nonuniform but 
static electric field because the electrons moving trans- 
verse to the electric field will not feel any change in sign 
of the electric field. 



V. DISCUSSION AND SUMMARY 

Spin accumulation observed by Sih et al. in two di- 
mensional electron gas corresponds to the value of 8 ~ 
0.1. This is in the large disorder limit. Although the 
applied electric field is uniform, the electronic inhomo- 
geneity in the system due to disorder may cause spatial 




FIG. 6: (Color online) In a spin-Hall bar geometry, an electric 
current j% is driven through a 2DEG contacted with the met- 
alic leads connected to a voltage source. The electric field E x 
with transverse momentum q y induces a novel nonvanishing 
bulk spin-Hall current j* with spin pointing pependicular to 
the plane of 2DEG. 



variation of the electric field in the system. The varia- 
tion of the electric field describes the presence of modes 
q. If transverse q is closer to 2mA, the contribution of 
these modes to the spin Hall conductivity is not negligi- 
ble. The presence of anomalous spin Hall current is then 
certain although the magnitude may be small since S is 
small in this experiment 4 . This contribution is intrin- 
sic because the spin orbit interaction is not disordered. 
Consideration of extrinsic mechanisms along with the in- 
trinsic spin orbit interaction provides finite ,35 spin Hall 
conductivity in presence of uniform electric field as well. 
However, any quantitative comparison of the anomalous 
spin Hall conductivity presented in this paper with the 
experiment 4 or with the contribution arising from extrin- 
sic effect 34,35 or with the spin accumulation across the 
edges in a ballistic system 36 is beyond the scope of the 
present study. Nevertheless our theory may be tested by 
applying a spatially varying electric field with the varia- 
tion along its transverse direction. The geometry for an 
experimental proposal to test the novel mechanism for 
'anomalous' spin Hall current is depicted in Fig. 6. 

The other studies at finite q in this Rashba 2DEG are 
the induction of spin-density by electromagnetic wave 37 , 
the response of the in-plane polarization 38 to the trans- 
verse electric field in a pure system and determination of 
density-density correlation function 39 at all q. 

In a cubic Rashba model which is relevant for two di- 
mensional hole gas, the intrinsic SHC is nonzero 16 , but 
the conserved spin Hall conductivity vanishes 2 for short 
ranged impurity potential. Therefore it is indeed inter- 
esting to look if these systems also have anomalous spin 
Hall current 40 like we have described here. 

In summary, we have determined spin Hall conductiv- 
ity at finite frequency and finite transverse wavevector in 
a disordered two dimensional electron gas with Rashba 
spin orbit interaction. Interestingly at zero or small fre- 
quencies, we have found an anomalous spin Hall conduc- 
tivity which resonates when the wavelength of the spatial 
variation of the electric field matches with the length of 
spin precession. The mechanism responsible for this is 
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the change in the direction of spin precession for elec- ishnan and K. Ray for stimulating discussions and M. 
trons moving perpendicular to the electric field when the Plctyukhov for sharing his calculation before publication, 
sign of the electric field changes due to spatial variation. 
This is primarily due to the change in sign of the out-of- 
planc component of spin. 
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